Numerical simulation of copper ablation by ultrashort laser pulses 
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Using a modified self-consistent one-dimensional hydrodynamic lagrangian fiuid code, laser ab- 
lation of solid copper by ultrashort laser pulses in vacuum was simulated to study fundamental 
mechanisms and to provide a guide for drilling periodic microholes or microgratings on the metal 
surface. The simulated laser ablation threshold is a approximate constancy in femtosecond regime 
and increases as the square root of pulse duration in picosecond regime. The ablation depth as a 
function of pulse duration shows four different regimes and a minimum for a pulse duration of ~ 12 
ps for various laser fluences. The influence of laser-induced plasma shielding on ablation depth is 
also studied. 



I. INTRODUCTION 

Micromachining with ultrashort laser pulses has been a 
rapidly advancing area of ultrashort laser applications in 
recent twenty years. For its high precise and practically 
melting free processing of nearly all kind of materials, 
it has been applied to many fields such as pulsed laser 
deposition (PLD), nanoparticle manufacturing, surgery 
and chemical analysis. It utilizes the properties of fem- 
tosecond or picosecond laser pulses to achieve an un- 
precedented degree of control in sculpting the desired 
microstructures internal or external to the material with- 
out much collateral damage to the surroundings [1] . The 
quality of ablated holes and patterns produced with ul- 
trashort laser pulses is much better than the quality 
of structures produced with nanosecond or microsecond 
laser pulses. Another new application of ultrashort laser 
pulses was reported that femtosecond laser processing 
could create a variety of colors on metal surface by laser- 
induced periodic microholes and could turned aluminum- 
like metals into golden appearance [^. This metal col- 
orization techinique will ultimately allow people to con- 
trol the optional properties of metal surface. 

To control and optimize the process of laser microma- 
chining, the understanding of the fundamental mecha- 
nisms of ultrashort laser ablation is required. The cou- 
pling of the ultrashort laser pulse to the metals is a com- 
plex process. Several theoretical models [ll y, |3] have 
been proposed to simulate ultrashort laser ablation of 
metals. Some other models such as electrostatic ablation 
[5| and Coulomb explosion [6] were also proposed but 
these models are not supposed to apply for the metals 
due to their high electron conductivity and fast thermal 
diffusion. During ultrashort laser ablation of aluminum, 
as described in Ref. [7] , clear separation between the ab- 
lated matter and the unablated target has been seen to 
occur through spinodal decomposition involving thermo- 
dynamic instabilities near the critical point of aluminum. 
Fast heating and phase transition leading to phase explo- 
sion seems to be a more realistic mechanism in metal ab- 
lation by ultrashort laser pulses and hydrodynamic mod- 
els can provide a consistent description of the processes 
of laser ablation of the metals when the laser intensity is 
not too high but appropriate for micromaching use. 



In this paper, we simulated the copper ablation by fem- 
tosecond and picosecond laser pulses with a modified hy- 
drodynamic code MED103 [3, Q. The dependence of the 
copper laser ablation as a function of laser pulse parame- 
ters (such as laser fluence, pulse duration, the number of 
pulses and repetition rate) and other effects which were 
involved in the process of copper laser ablation. With 
the aim to improve our understanding of the physical 
processes involved in laser-driven ablation and determin- 
ing the optimum physical condition for manufacturing re- 
quired surface structures, we have undertaken the devel- 
opment of a computer code treating both the processes of 
laser-driven ablation and subsequent plasma expansion in 
vacuum. A one-dimensional fluid model was employed in 
the computer code which include a self-consistent treat- 
ment of hydrodynamics, laser energy absorption, ther- 
mal conduction, equilibrium between electrons and ions 
temperature, fluid motion, shock waves generation and 
etc. Similar codes [3, [Ifl, [ll| have been developed to 
simulate ultrashort laser ablation, laser-driven thermonu- 
clear fusion, hot-electron generation and transport, laser- 
induced plasma spectroscopy in recent twenty years. A 
set of experimental conditions have to be carefully cho- 
sen in order to create required structures on the metal 
surface, so complete simulation is essential to the success 
of the experience. 

The following section describes the model used in the 
hydrodynamic code, and Sec. Ill presents the simulation 
results. A summary and discussion are given in Sec. IV. 



II. THEORETICAL MODEL 

Hydrodynamic codes are a useful tool for a self- 
consistent description of laser ablation processes when 
the laser intensity is not too high. At high laser intensi- 
ties collisionless and relativistic effects become dominant 
[3, 0, [l0l - |l2j . Particle in cell (PIC) simulations including 
collisions indicate the validity energy range of hydrody- 
namic simulations goes up to S'lA^ « 10^"^ {W/cm^)fi'm? 
where Sl is the laser pulse intensity and A is the laser 
wavelength [IQ] . 

The coupling of the ultrashort laser pulses to the metal 
target is a complicated process which is very different 



from those of semi-conductors or insulators because of its 
high thermal conductivities and its low melting temper- 
ature. The process involves the absorption of the laser 
pulse energy in the cold metal surface, the generation 
of high temperature and hydrodynamic pressure gradi- 
ents, the propagation of electron heat wave propagating 
into the bulk of the metal and the laser-induced plasma 
expansion etc. For femtosecond laser pulses, laser en- 
ergy is very rapidly absorbed by the electrons meanwhile 
lattices remain in original temperature. In metals, the 
conduction band is partially occupied by electrons. The 
electrons gain energy by inverse-Bremsstrahlung absorp- 
tion and nonlinear absorptions like photoionization and 
avalanche ionization to be excited to a higher energy level 
in the conduction band. At the end of the laser pulse 
the energy distribution of electrons may be noneqilib- 
rium and this nonequilibrium energy distribution ralax 
to a Fermi-Dirac distribution in sub-picosecond which 
can be properly characterized by a macroscopic quan- 
tity, i.e, temperature T. For picosecond laser pulses, the 
electron energy distribution generally reaches an equilib- 
rium state at the end of the laser pulses. During the 
laser ablation thermal process does not have much con- 
tribution because the duration of the whole process is 
shorter than the electron-lattice relaxation time which 
is typically several picoseconds. As a consequence of 
rapid laser energy deposition, irradiating metals with ul- 
trashort laser pulses results in strong noneqilibrium con- 
ditions between the electrons and the lattice. In this 
case electrons and lattice subsystem can be described by 
the classical two-temperature model (TTM) proposed by 
Anisimov et al. [l3|. The huge temperature gradients 
from nonequilibrium hot electrons to relatively cold lat- 
tices in a constant volume lead to the buildup of high 
compressive stresses which result in the development of 
a tensile component of the pressure wave that propagates 
deeper into the bulk of the target. The tensile stress in- 
creases with the depth until it reach the maximum value. 
If this stress could eventually overcome the mechanical 
strength of the target material, it could lead to the me- 
chanical separation and ejection of material, i.e., ablation 
occurs. Electrons and ions temperature eventually reach 
an equilibrium and then thermal expansion occurs. The 
metal surface layer on the order of the optical penetra- 
tion depth is ablated by electron emission, sublimation 
and transition to the plasma state. The remaining heat 
diffuses into the metal and leads to the emission of parti- 
cles and droplets due to the thermal boiling process. The 
material ejected could carry away most of the deposited 
laser energy, especially when working very close to the 
ablation threshold. Because of minimizing the amount of 
thermal energy diffused into the metal, ultrashort laser 
ablation greatly reduced the collateral damage to mate- 
rial surrounding the ablation zone. 

The code used here is a hydrodynamic Lagrangian fluid 
code, which includes ionization, collisional and resonance 
absorption, thermal transport, equation of state (EOS), 
hydrodynamics and a number of optional physics such 



as soft X-ray emission and lasing [8|, |9|. In the hydro- 
dynamic model, the fluid equations for the conservation 
of mass, momentum, and internal energy are solved with 
finite difference method. This model treats the electrons 
and ions as different subsystem characterized by internal 
energy U, temperature T, pressure p, specific heat ratio 
7 and so on. They are coupled together by the energy ex- 
change due to the electron-phonon or the electron-atom 
collisions. In the Lagrangian numerical scheme, the mass 
of each mesh is constant in time. The velocity u which 
defines the motion of the Lagrangian coordinates is ob- 
tained by solving the Navier-Stokes equation. The tem- 
perature is obtained by solving the energy conservation 
equation: 
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where U is the internal energies per unit mass and source 
term S is the rate per unit mass at which energy enters 
each subsystem. The source term Si and Se for ions and 
electrons respectively are wirtten as S'i — Hi — K + Q 
and Se — H,, + K + J + X where H represents the flow 
of heat due to thermal conduction, K = V ■ kVT / p is 
the rate of energy exchange between ions and electrons, 
J is the rate of Bremsstrahlung emission, X is t he rate 
of absorption of laser pules and Q is the rate of viscous 
shock heating. The thermal conductivity k is obtained 
from the model of L.Spitzer with correction factors. For 
the hot plasma we use Spitzer's formula for the collision 
frequency [14| : 
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where Z* is the ionization degree and ln{K) is the 
Coulomb logarithm. The non-LTE (non-local temper- 
ature equilibrium) time-dependent average model was 
used to calculate the electron energy obtained in the 
atomic processes of excitation, ionization etc. With 
this model, the rapidly heated plasmas induced by 
the ultrashort laser pulse is simulated. The elec- 
tron thermal heat conductivity KSpUzer is defined by 
K*nekBTe/ inieUJspitzer), whcrc K* is given by k* = 
(5e320/(37r) where 5 and e are the correction factors, 
which are weakly dependent on Z* and tabulated jl4|. 
The exchange rate K of energy between ions and elec- 
trons is given hy K = uneksiTi — Te)/p, where i/ 
is the exchange rate expressed as 3mekBijJspitzer/'mi- 
Laser absorption is assumed to occur via inverse- 
Bremsstrahlung at densities below the critical densities 
Pc = SaMraHmeUj\/ {Ze^) kg/m? at which the laser- 
induced plasma frequency equals the frequency ujl of 
the laser pulse. The resonance absorption is simply 
considered that twenty-five percent of laser energy hav- 
ing reached critical density is deposited in the plasma 
by resonance absorption. The distribution of the fast 
electrons produced by resonance absorption is approx- 
imately a Maxwellian distribution at temperature Th- 



The compilation of Th scales as (<^X^Y/^ for ($A^) > 
IQis Wcm-^^im^ and as [^X^f'^ for (^A^) < 10^^ 
Wcm~^IJLm^ where $ is the laser intensity arriving at 
the critical density surface. The laser pulses energy ab- 
sorption term X{r,t) is given by X{r,t) = PL{r,t)/dAI, 
where Pl {r, t) is obtained from laer energy absorption 
via inverse-Bremsstrahlung at densities below the critical 
density pc and resonance absorption. The EOS used in 
the code assumes the ions behaved as a non-degenerate 
perfect gas and the electrons behaved as a perfect gas 
which might either be non-degenerate, partially degener- 
are or fully degenerate. The shock wave created by the 
expanding plasma is consistently treated as for the tar- 
get. This model also includes a high- field corrections to 
inverse-Bremsstrahlung absorption, tunnel ionization, a 
non-local electron heat conduction and the modelling of 
coUisional in addition to recombination X-ray lasers. 



III. RESULTS 



to the fact that thermal penetration depth I = Dra is 
smaller than optical penetration depth 5 where D is the 
thermal conductivity and Ta is the ablation time. In this 
case threshold fluence Fth could be written as Fth — P^S 
where p is the density, J7 is the specific heat of evapora- 
tion. The optical penetration depth S can be expressed 
as i5 = c/wfc, where k is the imaginary part of the re- 
fractive index. The optical penetration depth does not 
rely on pulse duration, so the threshold fluences are ap- 
proximately the same for different pulse durations. The 
behavior of threshold fluence varying as the square root 
of the laser pulse duration r^ larger than 1 ps results from 
the threshold for ablation being controlled primarily by 
the thermal diffusion of the incident laser pulse enenrgy 
[J, [iJl- In this situation the heat penetration depth is 
much larger than optical penetration depth and thresh- 
old fluence can be written as Fth — P^l — P^D^a 
which is in direct proportion to Ta^l"^ . Note that in this 
case, tl ^ r^, the ablation time Tq « TeT^/ije + r^), 
where Tg and r^ are the electron cooling and the lattice 



Ablation threshold 



heating times, so threshold fluence varies as r^ 
is understandable. 
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The most popular industrial application of the ultra- 
short laser is micromaching such as drilling or cutting. 
Such a process requires pulse energy enough to give us 
laser fluence above ablation threshold. However, an in- 
teresting application called laser induced periodic sur- 
face structure (LIPSS) [l^ applied for colorizing metals 
happens near the ablation threshold while even below 
the ablation threshold the ultrashort laser pulses can in- 
duce amorphization on the target without ablation, so 
the precise value of laser ablation threshold is needed. 
Ultrashort laser ablation occurs when the laser fluence 
exceeds a certain threshold fluence which relies on the 
material properties, laser wavelength, pulse intensity and 
duration. The target is assumed to be ablated when a 
rapid decrease in the copper density profile appears. The 
ablation threshold was determined by the laser fluence 
when the copper surface density is changed rapidly and 
the simulated ablation depth exceeded 10 nm which was 
estimated to be the minimum ablation depth. Fig. 1 
shows the simulation results of the threshold fluence of 
ablation in copper sample by ultrashort laser pulse as a 
function of the pulse duration. The threshold fluence in- 
creases as laser pulse duration and shows two different 
regimes with a transition occuring between 1 ps and 5 
ps. For femtosecond pulse, the threshold fluence is ap- 
proximately a constant value of 0.237 J /cm?, while for 
pulse longer than ^ b ps the threshold fluence could be 

1/2 

represented as 0.184r^ . The simulated threshold flu- 
ence is in qualitative agreement with the simulation re- 
sults of Ref. jlj, lis] and the experimental data obtained 
for gold [J|, fused silica and calcium fluoride [l3|- The 
constancy of the threshold fluence for femtosecond laser 
pulse is in good agreement with the measurements made 
for copper sample [131 ■ The constant value of the abla- 
tion threshold for femtosecond laser pulse is basically due 
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FIG. 1. Threshold fluence vs pulse durations for copper ab- 
lation in vacuum (800 nm wavelength, 33/s — 100 ps pulse 
duration range). 



B. Ablation depth 

From theoretical analysis derived from one- 
dimensional two-temperature model [13| . the ablation 
depths can be expressed as 



L ^ Sln{F/Fi) {6 » 0- 



L ~ IMF/Fl^) (S « /). 



(3) 



(4) 



Fig. 2 shows the ablation depth of copper as a func- 
tion of the laser fluences obtained by our simulation with 



fixed pulse duration tl — 150 fs and laser wavelength 
A = 780 nm. The ablation depth generally shows two dif- 
ferent logarithmic dependence on the laser fluence which 
is in qualitative agreement with the reported experimen- 
tal results [1]. For laser fluence F less than 0.9 J /cm?, a 
good fit to the simulation results with equation (3) gives 
the threshold fluence Ff^^ ~ 31 mJ /cm? and optical pen- 
etration depth 5 ~ 10.4 nm which is in good agreement 
with standard value (5 = 13 nm, [19j . In this regime, 
the laser energy absorption coefficient is so low that not 
much energy transferred deeper than the skin depth of 
the copper target. This ablation regime is dominated 
by the optical penetration and the influence of thermal 
conduction can be negligible. One observes a transition 
occurring between 0.9 J /cm? and 2.0 J /cm?. In this 
regime, the influence of electronic heat diffusion becomes 
more obvious with the increase of laser fluence. At higher 
fluence {F > 2.0 J/cm?), the simulation results fitted 
with equation (4) gives the threshold fluence i^^^ ~ 1.470 
J/cm? and thermal penetration depth I ~ 126.8 nm. 
The obtained thermalpenetration depth is bigger than 
experimental results [l| because in that experiment the 
laser pulses are influenced by the walls induced by former 
pulses which decreased the effective penetration depth. 



duration reduces the ablation depth for the same laser 
fluence. In case of pulse duration longer than 10 ps, no 
logarithmic dependence on the laser fluence can be ob- 
served in the region F > 2 J /cm? . As described in Ref. 
ill], this effect can be explained by the hydrodynamic 
plasma expansion during the laser pulse, plasma shield- 
ing of the laser radiation, and increased heat-conduction 
losses. 
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FIG. 2. Copper ablation depth versus laser fluence for 150fs 
laser pulse. The dash and the dash-dot lines represent the 



most compatible fits, where F-^f^ 
1.470 J/cm^. 



0.031 J/cm^ and F^^ 



In Fig. 3, simulation results of the ablation depth for 
different pulse durations from 500 f s to 20 ps are sum- 
marized. There have no much variations between the ab- 
lation depths in case of 500 fs and 1 ps where the laser 
energy deposition is determined by the optical penetra- 
tion depth as explained before. In case of pulse duration 
T = 3 ps and r = 5 ps, the ablation depth values are obvi- 
ously lower than the ones of 500 fs and 1 ps for same laser 
fluence in the regime -F > 1 J /cm? . In this regime, the 
logarithmic dependence on the laser fluence still works, 
but the ablation depths and the threshold fluences are 
now dependent on the pulse duration. Increasing pulse 



FIG. 3. Copper ablation depth versus laser fluence for dif- 
ferent pulse duration: 500 fs, 1 ps, 3 ps, 5 ps, 10 ps, 20 
ps. 

Fig. 4 presents the simulation results of the copper 
ablation depth versus the pulse duration r with differ- 
ent laser fluences. The ablation depth decrease rapidly 
in the region r < 150 /s. In this region resonance ab- 
sorption plays an important role in the process of laser 
energy absorption and produces hot-electrons with huge 
energy penetrating deeply into the bulk of the copper tar- 
get. In this simulation situation the temperature of hot- 
electron is given by Th ~ (^A^/t)^/'^. With the same 
laser fluence, the bigger pulse duration r is, the smaller 
the hot-electron's temperature is which is equivalent to 
the smaller ablation depth. No significant change in abla- 
tion depth is observed in the range of laser pulse duration 
from 0.15 ps to 2 ps. The relative constancy in this region 
can be explained by the arguments used in the ablation 
threshold subsection. In the region 2 ps < r < 50 ps, the 
ablation depth decreases with the increase of the pulse 
duration. In this ps regime, the laser energy absorption 
is dominated by inverse-Bremsstrahlung, which results in 
a decrease in the amount of energy reaching the copper 
surface, so that the ablation depth decreases as well. One 
observes a minimum in the ablation depth occurs in re- 
gion 10 < r < 100 ps. In the region 50 ps < t < 100 ps 
the ablation depth increases a little up to 70 nm at first 
and then reaches a relative constant value. An analytical 
investigations in Ref. [20| indicated that the laser absorp- 
tion increased with A/A where A = \dzln{ne)\^ ^„ is the 
electron density gradient length. The ablation depth is 
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FIG. 4. Ablation depth as a function of the pulse duration 
for different laser fiuences of 5.4, 10.3, and 20.6 J/cm^ . 



proportion to the total absorbed laser energy in case of 
Gaussian pulses, so the increase of the ablation depth in 
this region can be explained by the the increase of the 
energy absorption in the laser-induced electron plasma. 



C. Multiplepulse 

Laser ablation efficiency can be optimized by using a 
prepulse to create a plasma profile on the surface of the 
metal before the main pulse [3| . This effect is determined 
by the timing between the prepulse and the main pulse 
because if the temporal separation is too long, laser en- 
ergy is absorbed in the plasma far away from the solid 
surface and does not contribute to the laser ablation. In 
this situation, ablation depth largely depends on the lag 
between the prepulse and the main pulse. The simula- 
tion results of ablation depth as a function of the lag 
between the prepulse and the main pulse are shown in 
Fig. 5. The general shape of the curves is in qualita- 
tive agreement with the experimental results described 
in Ref . [2l[ . Three different regimes of the laser ablation 
denoted with Latin alphabets depending on the lag be- 
tween the prepulse and the main pulse. The constancy of 
the ablation depth in region A where the lag is less than 
1.3 ps indicates that the double pulses are close enough to 
each other which could be considered as one pulse. With 
the lag longer than 8 ps, in region C, the main pulse is 
completely shielded by the prepulse-induced plasma. The 
region B of 1.3 ps ~ 8 ps corresponds to the transition 
regime with partial plasma shielding. 

We utilize upto five laser pulses to simulate the de- 
pendence of the ablation depth on laser pulse number in 
order to investigate the influence of plasma shielding on 
copper ablation in case of several laser pulses. Fig. 6 in- 
dicates the total ablation depth increase with increasing 
number of laser pulses, while ablation depth per pulse 



FIG. 5. Copper ablation depth as a function of the lag be- 
tween main pulse and prepulse. The used prepulse and main 
pulse parameters: pulse duration r = 150 /s, laser wave- 
length A — 1.0/i m. The fluence of prepulse is about 1.64 
J/crn^ and total laser fluence is about 18 J /cm?. 



decrease with increasing number of pulses. The delay 
between two pulses was set as to be near 3.0 ps and 
the fluence of each pulse was set to be near 8 J jcvr? . 
Experimental results presented in Ref. [24I indicate an 
approximate linear increase of ablation depth with in- 
creasing number of laser pulses. Our simulation results 
deviate from the linear dependence apparently. The rea- 
son is that the lag between pulses in our simulation is so 
short that the influence of plasma shieding is important, 
while in the experiment the lag between pulses was so 
long that the plasma shielding effect can be neglect and 
preceding pulses almost have no influence on subsequent 
pulses. Ultrashort laser drilled holes typically exhibit a 
conical shape due to the Gaussian intensity profile of the 
laser beam and enhanced thermal losses at the walls of 
the hole. As a consequence, the laser pulses are inter- 
acting with the walls at grazing incidence increasing the 
interaction area that decrease the effective fluence. 



IV. SUMMARY 

This paper presents the simulation results and physi- 
cal interpretations of copper ablation by ultrashort laser. 
The ablation threshold shows a approximate constant 

value for pulse duration less than 1 ps and increases as 

1/2 
0.184r^ for pulse duration tl > 2 ps after a transition 

occuring between 1 ps and 2 ps. This dependency re- 
lationship can be analytically explained by different ab- 
lation processes in the cases of pulse duration smaller 
or greater than the electron-lattice relaxaion time. In 
the former case, deposited laser energy during the pulse 
duration has no much difference compared to laser en- 
ergy leading to ablation after the complete of equilibrium 
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tal results reported in Ref. [1|. The dependence of ab- 
lation depth on pulse duration for a certain laser fluence 
clearly shows the efficiency of laser ablation of copper 
in femtosecond regime is higher compared to the one in 
picosecond regime. A experimental investigation of Ref. 
|23| indicate that laser ablation efhciency can be opti- 
mized by using a prepulse to create a preplasma for the 
main pulse, but our simulation of copper laser ablation 
shows a continuing reduction of ablation depth as the in- 
crease of the lag between the prepulse and main pulse. 
Because of the plasma shielding effect, the laser energy 
of subsequent pulses can't completely contribute to the 
laser ablation of copper if the lag between two pulses is 
too short. In general we should set the lag between two 
pulses larger than dozens of nanoseconds when the effect 
of plasma shielding has no much influence. 



FIG. 6. Total ablation depth and ablation depth per pulse as 
a function of the number of pulses in different pulse duration 
100 fs, 150 fs, 500 fs. 
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